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1 Introduction 

During the last decade the considerable progress in the rigorous derivation in scaling limits 
of quantum kinetic equations, in particular the nonlinear Schrodinger equation [1]- [16] and 
the Gross-Pitaevskii equation [H]- [19] as well as the quantum Boltzmann equation [20] is 
observed. 

In paper f2T] we established that, if initial data is completely defined by a one-particle 
marginal density operator, then all possible states of infinite-particle systems at arbitrary mo- 
ment of time can be described by the generalized quantum kinetic equation within the frame- 
work of a one-particle density operator without any approximations. The aim of this paper is 
to construct the mean field (self-consistent field) asymptotics of a solution of the initial-value 
problem of the generalized quantum kinetic equation and to extend this result on case of the 
kinetic evolution in the presence of initial correlations of quantum particles. 

We shortly outline the structure of the paper and the main results. At first in Section 

2 we formulate some definitions and preliminary facts on the description of quantum kinetic 
evolution. Then the main results on the mean field scaling limit of a solution of initial-value 
problem of the generalized quantum kinetic equation and the marginal functionals of the state 
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are stated. In Section 3 we prove the main results. We established that the constructed asymp- 
totics of a solution of the generalized quantum kinetic equation is governed by the quantum 
Vlasov kinetic equation for the limit states and the limit marginal functionals of the state are 
products of a solution of the derived Vlasov kinetic equation. In Section 4 we consider some 
consequences and generalizations of the obtained results. In particular, we extend the quan- 
tum kinetic equations on case of the evolution of particle states in the presence of correlations 
at initial time. Finally in Section 5 we conclude with some observations and perspectives for 
future research. 

2 Mean field dynamics of quantum many-particle sys- 
tems 

We adduce some definitions and preliminary facts about the description of quantum dynamics 
within the framework of a one-particle density operator governed by the generalized quantum 
kinetic equation. Then the main results about its mean field scaling limit is formulated in case 
of quantum particles obeying the Maxwell- Boltzmann statistics. 

2.1 The generalized quantum kinetic equation 

We consider a quantum system of a non-fixed (i.e. arbitrary but finite) number of identical 
(spinless) particles obeying Maxwell-Boltzmann statistics in the space M.'^. We will use units 
where h = 2nh = 1 is a Planck constant, and m = 1 is the mass of particles. 

Let T-Lhe a. one-particle Hilbert space, then the n-particle spaces Hn, are tensor products of n 
Hilbert spaces H. We adopt the usual convention that H'^^ — C. We denote by — 0^o^" 
the Fock space over the Hilbert space Ti. 

The Hamiltonian of n-particle system is a self-adjoint operator with domain T>[Hn) C 1-Ln- 

n n 

i7„ = J]K(i)+6 J2 *(^l'^2), (1) 

i=l ii<i2=l 

where K{i) is the operator of a kinetic energy of the i particle, 12) is the operator of a two- 
body interaction potential and e > is a scaling parameter. The operator K{i) acts on functions 
ipn, that belong to the subspace Ll{W"-) C T>{Hn) C L'^{W"') of infinitely differentiable func- 
tions with compact supports according to the formula: K{i)ipn — —^Ag^ipn- Correspondingly 
we have: $(^l,^2)'0n = ^(9«i)fe)'0n, and we assume that the function $(gji,gj2) is symmetric 
with respect to permutations of its arguments, translation-invariant and bounded function. 

Let = 0^o'^^('^") space of sequences / = (/o, /i, ...,/„,.. .) of trace 

class operators /„ = . . . ,n) G and /o G C, that satisfy the symmetry condition: 

/„(!, . . . ,n) — fn{ii, ■ ■ ■ ,in) for arbitrary {ii, . . . , i„) G (1, . . . , n), equipped with the norm 

00 00 

n=0 n=0 



Mean field asymptotics 



4 



where Tri_...^„ are partial traces over 1, . . . ,n particles. We denote by 2l{J-'-u) = ©^o-^ol^n) 
the everywhere dense set of finite sequences of degenerate operators with infinitely different iable 
kernels with compact supports [22] . 

On the space 2}{T-h) we define the group G{—t) = ®'^=QGn{—t) of operators of the von 
Neumann equations 

6^n(-t)/n = e-**^"/„e**^". (2) 

On the space 2^{J-'-u) the mapping ([2]): t Q{—t)f is an isometric strongly continuous group 
which preserves positivity and self-adjointness of operators. For G £o('H„) there exists a 
limit in the sense of the norm convergence on space £}{'Hs) by which the infinitesimal generator 
of the group of evolution operators ([2]) is determined as follows 

lim - /„) = -i{HrJn - fnHn) = -J\fnfn, (3) 

where if„ is the Hamiltonian ([1]) and the operator: —i{Hnfn — fnHn) is defined on the domain 
V{Hn) C Tin- We denote by {—Afmtih j)) the operator 

(-A/'int(^, j))/n = fn " j)) (4) 

defined on the subspace £q 

Let us denote Y = {1, . . . , s), X \ Y = (s + 1, . . . , s + n) and {¥} is the set consisting of one 
element Y = {1, . . . ,s), the mapping 9 is the declusterization mapping defined by the formula: 
9{{Y}, X\Y) = X. We define the (1 + ri)t/i-order {n > 0) cumulant of groups of operators ([2]) 
as follows [23] 

Qi,^n{t,{Y},X\Y)= Yl (-I)I^I-^(IPI-I)! nfe)lH'^TO)' (5) 

P :{{¥}, X\Y)=[J^X, X,CP 

where is the sum over all possible partitions P of the set {{Y} , X\Y) = {{Y}, . . . , s+n) 
into |P| nonempty mutually disjoint subsets Xi C {{Y},X \ Y), for example, 

^,{t,{Y}) = g,{-t), 

2l2(t, {Y},s + 1) = Gs+ii-t, Y,s + 1)- g,{-t, Y)gii-t, s + l). 

We indicate some properties of operators ([5]). If n = 0, the generator of the first-order 
cumulant 2ti(f, {Y}) = Qs{—i) for fs G 2l{'Hs) C 2^{'Hs), in the sense of the norm convergence 
on space £^(7/3), is given by the operator 

limJ-{Ql,it,{Y}) - I)f, = -KU 

In the case n = 1 we have in the sense of the norm convergence in £^(7/^+1) 

1 " 
lim - 2t2(t, {Y}, s + = -eJ2 -Mntl^, s + 

i=l 
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where the operator {—Mmtih s + 1)) is defined by formula (jlj), and for n > 1 in the sense of 
the norm convergence in ii^('Hs+n), as a consequence that we consider a system of particles 
interacting by a two-body potential, it holds 

lim^2li+„(t)/^+„ = 0. 

We consider the mean field (self-consistent field) asymptotic behavior of a solution of the 
Cauchy problem of the generalized quantum kinetic equation |f2Tj 



^Fi(t,l) = -A/-i(l)Fi(t,l)+ (6) 

oo ^ n+2 

+eTr2 ( - Mntil, 2)) Tr3,...,„+22Ji+n (t, {1, 2}, 3, . . . , n + 2) J] F,{t, t), 

n=0 ' 1=1 

Fi(t,l)|,=o = F°(l). (7) 

In kinetic equation ([6]) the {n + l)t/i-order generated evolution operator QJi_|.„(t), n > 0, is 
defined as follows (in case of {Y} = {1, 2}) 

fc=0 n.i=l nfc=l ^ '^^ 

x2ii+„_„,_..._nfe(t, {y}, s + 1, . . . , s + n - ?2i - . . . - rifc) x 

fc ^ s+n-7ii-...-nj ^ 

xH E E n W,i(M.,,^z,), 

|Dj I < s + n — rii — ■ ■ • — rtj 



where ^Dj Zj=u, ^i ^^'^ possible dissections of the linearly ordered set Zj = 

[s + n — rii — . . . — rij + 1, . . . , s + n — rii — . . . — Uj^i) on no more than s + n — rii — . . . — rij 
linearly ordered subsets, and we denote by 2li+„(f) the (1 + n)-order cumulant of the groups of 
scattering operators 

n 

Gn{t)=Qn{-tA,...,n)Wgi{t,i), n>l. (9) 

i=l 

For example, 

QJi(t,{r})=2ii(t,{y}), 

s 

Qj2(t, {1^}, s + 1) = 2i2(t, {y}, s + 1) - 2ii(t, {y}) 2t2(t, z, s + 1). 



i=l 



The collision integral series in kinetic equation (|6]) converges under the condition that: 
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The global in time solution of initial- value problem ©-(IT]) is determined by the following 
expansion [21] 

oo ^ n+1 

Fi(t, 1) = 5^ Tr2,...,i+„ 2ti+„(t, 1, . . . , n + 1) J] F°(z), (10) 

n=0 ■ j=l 

where 2ti+n(t) is the (1 + n)th-OYdeT cumulant ([5]) of groups of operators ([2]). The series (fTOil 
converges under the condition that: ||F°||£i(^) < e~^°(l + e~^)~^. 

In case of initial data given in means of a one-particle density operator the evolution of 
all possible states of quantum many-particle systems is described by a solution of initial-value 
problem of the generahzed quantum kinetic equation ([n])-(IZl) and a sequence of explicitly defined 
functionals of a solution of this generalized kinetic equation 

oo ^ s+n 

F,{t,Y I Fi(t)) = ^-Tr,+i,...,,+„Q3i+„(t,{r},X\y) J]Fi(t,z), s > 2, (11) 

n=0 ■ i=l 

where the {n + l)t/i-order {n > 0) generated evolution operator QJi+„(t) is defined by formula 
dH]). Marginal functional series ( ITTl) converges under the condition that: ||Fi(t)||£i(^) < e"*^^"*"^^^. 



2.2 The mean field limit theorems 

The mean field scaling limit of a solution of the initial-value problem of generalized kinetic 
equation ([6]) is described by the following limit theorem. 

Theorem 1. Let there exists the limit /° G -2^(7/) of initial data ^ 



then for finite time interval t G (— to,^o)) where to = (2 ||$||£(«2)ll/i IUh^)) ^' ^here exists the 
following limit of solution ^^) of the generalized quantum kinetic equation ^ 

lim||eFi(t)-/i(t)||^,(^j = 0, (12) 
where the limit one-particle marginal operator fi{t) is represented in the form 

OO „ „ 

/l(t,l) = ^ / dtl... / dt„Tl2,...,l+nGl{-t + tul){-Mnt{l,2)) X (13) 
n=o{ { 
2 n 

X JJ^l(-tl+t2,il)--- n +^n>in-l) X 

ii=i in-l=l 
n 1+n 1+n 



$^(-Mnt(^n,l+n)) llgi{-tn,Jn)llf'l 



°^Z). 
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For bounded interaction potentials series ( fT3l) is norm convergent on the space 
under the condition 

and for initial data /{* G ilJ('Hi) the scaling limit operator is a strong solution of Cauchy 
problem of the quantum Vlasov equation 

^A(t,l) = -M(l)/i(t,l) + Tr2(-Ari,,(l,2))/i(t,l)A(t,2), (14) 

/i(t)l*=o = /?. (15) 

Since a solution of initial- value problem ([n])-(IZl) of the generalized kinetic equation converges 
to a solution of initial- value problem (fT^ -( IT5|) of the quantum Vlasov kinetic equation as (fT2|) . 
for marginal functionals ( ITTj) we establish 

Theorem 2. Under the conditions of Theorem 1 for functionals 07]) it holds 

s 

lim \\e'F,{t, 1, . . . , s I Fi(t)) - J] /i(t, j) Li(«^) = 0, 

i=i 

where the operator fi{t) is defined by series l[T^) . 

This statement means that in the mean field scaling limit a chaos property preserves in time. 
In section 4 these theorems are extended on the case of kinetic evolution of quantum states 
in the presence of correlations of particles at initial time. 

3 The mean field limit of a solution of the generalized 
quantum kinetic equation 

We construct the mean-field scaling limit of a solution of initial-value problem of generalized 
kinetic equation (|6]) and marginal functionals of the state (ITTi) and prove stated above limit 
theorems. On the basis of obtained results we consider the problem of the justification of the 
nonlinear Schrodinger equation. 

3.1 Preliminaries: cumulants of asymptotically perturbed groups of 
operators 

For asymptotically perturbed first-order cumulant ([5]) the following statement is true [22] . 

Lemma 1. If fg G S^^iJ-Ls), then for arbitrary finite time interval for the strongly continuous 
group ^ it holds 

s 

lim||^.(-t)/.-n^iH>^U|U«,) = 0, 
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Proof. If an interaction potential is a bounded operator, then for fs G £^(7/5) an analog of the 
Duhamel formula for group ([2]) holds 

s 

{gs{-tA,...,s)-l[g,{-t,i))f,= 
1=1 

t 

/s s 
drl[gi{-t + r,l){- J2 Mnt{i,j))gs{-r)fs. (16) 
1=1 i<i=i 

Indeed, the Duamel equation ( !T6|) is valid for fs G illiT-Ls) C ii^('Hs). Since the operators from 
both sides of this equality are bounded and the set £,l{Hs) is everywhere dense set in the space 
£^('Hs), equality ( ITB]) holds for arbitrary fs G ii^('Hs). We note that the integral in (ITB]) exists 
in strong sense and the operator YU^iQii-t + r, /)( - Z]i<j=i-Mnt(«,i))^s(-T)/s is strongly 
continuous over r for every fs G £^(7{s), and hence it is integrable. 
Therefore the validity of lemma follows from the estimate 

s 

II {Qsi-t, 1, . . . , 5) - n Gii-t, o)/.Li(«,) < 
1=1 

t 

/s s 
dT\\l[g,{-t + T,l){- ■^int(^.j'))6^s(-r)/.||^,(^^) < 

< etsis-l)\m^^n2)\\fshHns)- 



□ 

In consequence of Lemma 1 for first-order cumulant ([5]) of scattering operators ([9]) the 
equality holds 

lim||^,(t)/, - /s||£i(^^) = 0. 



e-s-O 



Correspondingly, for cumulants ([5]) of asymptotically perturbed groups of operators we have: 

Lemma 2. If f s+n G 2}{'Hs+n), then for arbitrary finite time interval for the (1 + n)th-order 
cumulant of strongly continuous groups ^ it holds 

lim||li-2li+„(t,{y},X\F)/,+„- 



dti... / ^i^„JJ^;l(-^ + ^l,J)5Z(-A/'i„t(^l,s + l)) ^^l(-^l + ^2,Jl)••• 
s+rt— 1 s+n— 1 s+n 
■ W Gl{-tn-l + tn, jn-l) ( " -^mt (^n , S + n)) JJ ^1 , j„) 1 1 ^^^^^^ = 0. 
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The validity of this lemma follows from an analog of the Duhamel formula for cumulants ([5]) 
of strongly continuous groups ([2]) 

— rTfs+l,...,s+n2ll+n(^, {y}, X \ Y)fs+n = 

n\ 

t tn-1 ^ 

= j dh... j dtnTrs+l,...,s+nGsi-t + ^l) J] ( " Mntih, S + 1)) X 

*i=^ 

s+n— 1 

xGs+l{-ti +t2) . . . Gs+n-li-tn-1 + tn) ^ ( " Mnt{in, S + n))Q s+n{-tn) f s+n- 

and it is proving similar to previous lemma. 

We give instances of analogs of the Duhamel equation for cumulants (|5]) of scattering opera- 
tors ([9]) and generated evolution operators ([8]). If an interaction potential is a bounded operator 
and fs+i G £}{l-Ls+i)^ then for the second-order cumulant 2l2(t, {^}, s -f 1) of scattering opera- 
tors dl]) an analog of the Duhamel equation holds 

%2{t,{Y},S+l)fs+l= / dTgs{-T,Y)gi{-T,S + l)Y,{-J^U^US+l))x 

ii=i 

s+l 

X^.+l(r - t, y, S + 1) JJ Qi^T, l2)fs+l, 

12 = 1 

and, consequently, for the second-order generated evolution operator 2J2(^, s + l) we have: 

s 

2J2(t, {Y},s + = {Mi, 5 + 1) - {¥}) Mt, tus + = 

ii=i 

i-t s 

= / dTg,{-T,Y)g,{-T,s + i){y2i-^intiii,s + i))g,+iiT~t,Y,s + i)- 

s s+l 

-gs{T - t, Y) J2i-^int{^u s + l))^2(r - t, zi, s + 1)) n ^2)/s+l. 



il=l 12=1 



Then, according to Lemma 2, i.e. formulas of an asymptotic perturbation of cumulants of 
groups of operators, and definition (|8]) of the generated evolution operators, we establish 

limJ|(2Ji(t,{y})-/)/,||,,,_ = 0, 



and in the general case the following equalities hold 



Jim II 1 ^,^4t, {Y},X\ = 0, n>l. 
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3.2 The proof of the Hmit theorem 

We give a sketch of tlie prove of Tlieorem 1. 

In view that the series for eFi{t) converges, in the sense of the norm convergence on the 
space H^^H), under the condition that: 

t<tQ= (2||$||£(^,)||eF°||£i(«))~\ 

then for t < to the remainders of solution series ffTOj) and f|T3|) can be made arbitrary small 
for sufficient large n = uq independently of e. Then, according to Lemma 1, Lemma 2 and 
definition ([HD, for each integer n every term of these series converge term by term. 

Let us construct an evolution equation, which satisfies expression (fT3l) . We prove that it is 
a solution of initial- value problem (fT^ - (fT5|) of the quantum Vlasov kinetic equation. 

Taking into account the validity of equality ([3]), we differentiate expression (fT3l) over the 
time variable in the sense of pointwise convergence of the space £^(7^) 

|A(t,i) = -Ar(i)A(t,i)+ (17) 



t 

oo 



+Tr2(-A/'int(l,2))^ dh... / dtnTTs,...,n+2Y[Gl{-t + X 
n=0 Q Q ii = l 

2 3 n+1 

{-Mnt{h,3)) n^l(-^l+^2,il)... n Gl{-tn+tn,in) X 





2 3 n+1 

X 

kl = l jl=l jn=l 

n+1 n+2 n+2 

X Y,{-MUkn,n + 2)) n^^l(-^n,Jn)nA°W- 

kn=l jn = l i=l 



Using the product formula for the one-particle marginal density operator fi{t,i) defined by 
series (IT^ 



OO 



= ^ / (iti . . . / (i)f:„Trfc+i,...,fc+„ J]^ + ti 

j=l n=0 Q Q «i = l 

k fc+1 fc+n-1 

fci=i ii=i j„=i 

fe+n— 1 fc+n fc+n 



where the group property of one-parameter mapping ([2]) is applied, we express the second sum- 
mand in the right-hand side of equality (fTTl) in terms of operators ni=i/i('^; 0; ^^cl consequently, 
we derive kinetic equation ([1] 
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3.3 The limit marginal functionals of the state: the propagation of 
a chaos 

To give a sketch of the prove of Theorem 2 we represent marginal functionals of the state f ill I) 
in terms of the marginal correlation functionals Gs(t,Y \ s > 2, namely 

P:Y = \JiYi YiCP 

where ^p.y=u y. is the sum over all possible partitions P of the set y = (1, . . . , s) into |P| 

nonempty mutually disjoint subsets Yi C Y. The marginal correlation functionals Gs[t,Y \ 
Flit)), s > 2, are represented by the following expansions [23] 

oo _ s+n 

G,{t, Y I Fi(t)) = 5^ ^ Tr,+i,...,,+„QJi+„(t, e{{Y}),X \ Y) J] F,{t, z), s>2. (18) 

n=0 ' i=l 

In series f|T8|) it is introduced the notion of the declusterization mapping 6 : {Y} — )■ Y defined 
above. Hence in contrast to expansion f ITT]) the n term of expansions flTH]) of the marginal 
correlation functional Gs(t,Y \ Fi{t)) is governed by the (1 + n)t/i-order generated evolution 
operator ([8]) of the [s + n)t/i-order cumulants of the scattering operators. For example, the 
lower orders generated evolution operators 2Ji+„(t, 9({Y}), X \ Y) , n > 0, have the form 

QJi(t,^({F}))=2i,(t,^({y}), 

s 

^2{t, 0{{Y}),s + 1) = 2l,+i(t, 9{{Y}),s + 1) - 2i,(t, 9{{Y})) ^ m^it, i,s + l), 

1=1 

and in case of s = 2, we have 

QJi(t,^({l,2})) = ^2(t,l,2)-J. 

The statement of Theorem 2 is true in consequence of the validity for functionals f|T8|) the 
equalities 

lim||e^a(t,y|Fi(t))||^,(^^) = 0, s>2. 

In view of the structure of the generated evolution operators 2Ji+„ {t, 9{{Y}), X \ Y) , n > 0, of 
functionals of the state (ITSj) the last equality is true according to Lemma 1 and Lemma 2. 

3.4 Mean field quantum kinetic equations 

If we consider pure states, i.e. /i(t) = \i/jt){i't\ is a one-dimensional projector onto a unit vector 
\ipt) € H or in terms a kernel of the marginal operator fi{t): fi{t,q,q') = ip{t,q)ilj*{t,q'), then 
the quantum Vlasov kinetic equation reduces to the Hartree equation 

i-i:{t,q) = --Agij{t,q)+ / dq'^q - q')^, q')\'ij{t, q) . (19) 
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Moreover, for pure states, if it holds 

limlleFf- iV'o) (V'ol =0, 
the statement of Theorem 2 reads 

lim II e^F,{t I - jj^.^^^^ = 0, 

where l"^*) is the solution of the nonlinear Hartree equation (fT9|) for initial data {ipo)- 

We remark that in case of a system of particles, interacting by the potential which kernel 
is the Dirac measure = S{q), the Hartree equation fll9l) is reduced to the cubic nonlinear 
Schrodinger equation 

= -^A,i;{t,q) + \ij{t,q)\^^{t,q). 

The obtained results can be generalized on systems of quantum particles interacting via 
many-body potentials, i.e. systems with the Hamilton operators 

n n n 

1=1 k=2 ji<...<ifc=l 

where the operator is an operator of a fc-body interaction potential. 
In this case the generalized quantum kinetic equation (E]) has the form [2T] 

-F,{t, 1) = 1) + EE^'7;^^Tr,,...,„^,(-<t-^^))(l, . . . , (20) 

n=l k=l ^ 

n+1 

k + l)QJi+n_fe(t, {1, . . . , A; + 1}, A; + 2, . . . , n + 1) JJ Fi(t, z), 

i=l 

where the operator QJi+„_fc(t) is the (1 + n — A;)t/i-order generated evolution operator ([8]) and 

(-A/-(';))(i,...,fc)/. = -z(<i.(^)/.-/.<i.(^)). 

In case of a A;-body interaction potential the collision integral of the generalized quantum kinetic 
equation fl20|) is given by the norm convergent series under the condition that: < 

The mean field scaling limit fi{t) of solution ffTOj) of initial- value problem of the generalized 
kinetic equation f l20|) is described by the limit theorem similar to Theorem 2 and for initial 
data /° G -Co(^) ^ strong solution of the Cauchy problem of the following Vlasov quantum 
kinetic equation 

J oo ^ n+1 

-hit, 1) = -M{l)h{t, 1) + E ^ Tr2,...,„+i(-<r'^)(l, . . . , n + 1) n hit^ 0> 

n=l 1=1 
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Then for a many-body interaction potential the Hartree equation takes the form 



d 1 



n+l 



n=l ^' ^ i=2 

and correspondingly, we can derive the nonlinear Schrodinger equation with the 2n — 1 power 
nonlinear term. 

4 The kinetic evolution involving initial correlations 

One of the advantages of the developed approach is the possibility to construct the kinetic 
equations in scaling limits in the presence of correlations of particle states at initial time. 

We extend obtained results on case of quantum systems of particles which initial data spec- 
ified by initial correlations, for instance, correlations characterizing the condensate states of 
particles. 

4.1 Quantum kinetic equations in the presence of initial correlations 

We will consider initial state which is given by the following sequence of marginal density 
operators: 

2 n 

F(0) = (1, F°(l), <?2(1, 2) n ^i^^), • • • , 9n{l, . . . , n) n F^i^), ...), 

1=1 i=l 

where the bounded operators gn G n > 2, are specified initial correlations [24]- Such 

initial data is typical for the condensed states of quantum gases, for example, the equilibrium 
state of the Bose condensate satisfies the weakening of correlation condition with the correlations 
which characterize the condensed state [25] . 

In this case the one-particle density operator Fi{t) is governed by the following generalized 
quantum kinetic equation [26] 

|F,(t,l) = -Ar(l)Fi(t,l)+ (21) 



n=0 

?0/ 



oo ^ n+2 

+eTr2(-A/;nt(l, 2)) Tr3,...,„+2 C5i+„(t, {1, 2}, 3, . . . , n + 2) J] Fi{t, i), 



F,{t,l)\,=o = F^{l), (22) 
where the (1 + n)t/i-order generated evolution operator 0i+„(t), n > 0, is defined by the 
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following expansion: 

n n n-ni-...-nfc_i 

&,ut.{y}.x\Y) = n\Y^{-ifY.--- E 7 (23) 

^-^ ^-^ ^ — ^ [n — rii — . . . — nh]\ 

x2li+„_„i_..._„^(t, {F}, s + 1, . . . , s + n - ni - . . . - nfc) X 

k ^ s+71-ni-...-nj ^ 

Dj I < s + n — 711 — ■ ■ ■ — nj 

In formula (1231) we denote by XId z =u ^'^^^ possible dissections of the linearly 

ordered set Zj = {s + n — rii — . . . — rij + 1, . . . , s + n — rii — . . . — fij^i) on no more than 
s + n — rii — . . . — rij linearly ordered subsets and we introduce the (1 + n)th-oTdeT scattering 
cumulants 

s+n 

{Y}, x\Y) = 2ti+„(-t, {Y},x\ x \ r) J] 2ii(t, t). 



1=1 



For example, 

(3^{t,{Y})=^^{t,{Y}) = 

s 

= ^^{-t,{Y})g^{{Y})l[^^{t,^), 

1=1 

and 



(32it,{Y},s + l) 



s+l 



2t2(-t, {Y},s + l)92i{Y}, s + 1) n Mt, - 

i=l 

s s 

-M-t, {y})9ii{y}) n E ^2(-^, ^, s + 1)^2(^, s + i)2ii(t, ^)2ti(t, s + 1), 



j=i j=i 



where it is used notations accepted above. 

The global in time solution of initial- value problem fl2T|) - fl22|) is determined by the series [26] 



n+l 



F,{t, 1) = E ^ Tr2,...,i+„ 1, . . . , n + l)(7i+„(l, n + 1) J] (24) 

n=0 ' i=l 

where 2ti+„(t) is the (1 + n)th-oTdeY cumulant of groups of operators ([2]) and the operators 
9i+n, n > 0, are specified initial correlations. The series (^M converges under the condition 
that: <e-i0(l + e-9)-i. 

Correspondingly, the marginal functionals of the state are represented by the following ex- 
pansions: 

oo ^ s+n 

F,{t, Y I Flit)) = E Tr,+i,...,,+„ l5i+„(t, {Y},X\Y)l[ Fi(t, z), s > 2, (25) 

n=0 ■ i=l 
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where generated evolution operators of these functionals are defined by formula ( l23l) . 

Thus, the coefficients of generalized quantum kinetic equation ( 12T|) and generated evolution 
operators ( 123|) of marginal functionals of the state (125|1 are determined by the operators specified 
initial correlations. 

4.2 The mean field evolution of initial correlations 

In case of initial state involving correlations for generated evolution operator fl23l) of asymptot- 
ically perturbed groups of operators in the mean field limit the following equality is valid 

lim||l©i+„,(t,{y},X\y)/,+„|| =0, n>l, (26) 



and in case of first-order generated evolution operator (1231) we have, respectively 

s s 

lim II (0i(t, {Y}) - H g^{-t, ^^)g^i{Y}) J] Giit, ^2)) fs\\^.^^^^ = 0. (27) 



ii=i 



In view that under the condition that: t < to = (2 ||$||f;(-^2) ||e F^H^ic^)) ^, the series for 
eFi{t) is norm convergent, then for t < to the remainders of solution series (^M can be made 
arbitrary small for sufficient large n = no independently of e. Then, using stated above asymp- 
totic perturbation formulas, for each integer n every term of this series converges term by term 
to the limit operator fi{t) which is represented by the following series 

/i(t,l)= (28) 

CO „ „ z 

= dti... / dtnTl2,...,l+nQl{-t + ti,l){-UUl.2))\[Gi{-ti+t2,3l)... 

n=0 i { n=l 

n n 1+n 

... H g;i(-t„_i-i-t„,^„_i)^(-A/;,,(^„,i + n)) 

1 ( ) jn ) ^ 

in-l = l in = l jn = l 

1+n 

xg,^^{l,...,n+l)l[f,{i). 

1=1 



For bounded interaction potentials series (!28l) is norm convergent on the space £^('H) under 

the condition: t < to = {2 ||$||£(W2)||/?||£i(w))~ • 

Thus, if there exists the limit G -2^(7/) of initial data ([7j), namely 



then for finite time interval t G {—to, to), where to = (2 ||$||£(-K2)II/i IliiMW)) ^' there exists the 
mean field limit of solution expansion fl24|) of the generalized quantum kinetic equation fl2T]) : 



(29) 
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where the operator /i (t) is represented by series ( l28l) and it is a solution of the Cauchy problem 
of the modified Vlasov quantum kinetic equation 

lf,(t,l) = -Afil)Mt,l)+ (30) 

2 2 

+iY2(-Ari„t)(i, 2) H g^{-t, zi)(?i({i, 2}) n z2)/i(t, i)/i(t, 2), 

h = l 12=1 

/i(t)l*=o = /{'. (31) 

Since a solution of initial- value problem (I^T]) - fl22l) of the generalized kinetic equation con- 
verges to a solution of initial- value problem (1301) - (^1]) of the modified quantum Vlasov kinetic 
equation as fl29l) and equalities (|26|) and fl271) hold, for marginal functionals of the state fl25l) 
we establish 

lim||e^F,(t,F I Fi(t)) - J] g^{^-t,i,)g,{{Y}) J] 6^i(t, ^2) J] /i(^'^')Lmw.) = °- 



n=l 12=1 



This equality means the propagation of initial correlations in time in the mean field limit. 

Let us consider the pure states, i.e. the operator fi{t) = \'i/jt){4't\ is a one-dimensional 
projector onto a unit vector {ipt) G "H and its kernel has the following form: fi{t,q,q') = 
ijj{t,q)ilj*{t,q'). Then, we remark that in case of a system of particles, interacting by the 
potential which kernel $(g) = S{q) is the Dirac measure, the modified quantum Vlasov kinetic 
equation (1501) reduces to the Gross-Pitaevskii-type kinetic equation 

'TrMt^l) = -T^^gi^it^Q) + / dq'dq"b{t,q,q;q',q")i{j{t,q")r{t,q)ij{t,q), (32) 



ar^'^' 2 

where the coupling ratio b{t,q,q]q',q") of the collision integral is the kernel of the scattering 
length operator n?i=i Gii-t, «i)fei({l, 2}) n?2=i ^2)- 

Observing that on the macroscopic scale of the variation of variables, groups of operators 
(j2]) of finitely many particles depend on microscopic time variable £~^t, where e > is a scale 
parameter, the dimensionless marginal functionals of the state are represented in the form: 
Fs{e~^t I Fi{t)). As a result of the formal limit processing e — )■ in collision integral (!32|) . we 
establish the Markovian kinetic evolution with the corresponding coefficient b{e~^t). 



5 Conclusion 

In the mean field scaling limit we derived the quantum Vlasov kinetic equation and correspond- 
ingly, the Hartree equation (or the nonlinear Schrodinger equation) for pure states of quantum 
systems of particles obeying the Maxwell-Boltzmann statistics. In particular, in case of two- 
body interaction potentials it is the evolution equation with the cubic nonlinear term and in 
case of n-body interaction potentials the Hartree equation contains the 2n — l power nonlinear 
term. The obtained results can be extended to quantum systems of bosons and fermions. 

The mean field scaling asymptotics of a solution of the generalized quantum kinetic equation 
of many-particle systems in condensed states has been also constructed. We note that one more 
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approach to the construction of the kinetic equations in the mean field limit, in the case of the 
presence of correlations at initial time, can be developed on basis of the description of the 
kinetic evolution in terms of marginal observables [27]. 
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